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Experiments with imaginary chemical potential and Glasgow method carried out in two interrelated models 
- four dimensional QCD in the infinite coupling limit, and one dimensional QCD - support the point of view 
expressed by the title. 



The profound differences of dense, cold mat- 
ter, and hot, moderately dense, hadron gas might 
be reflected by different performances of numer- 
ical approaches to finite density QCD. In partic- 
ular, the enhancement of the fluctuations near 
Tc should give better chances either to the Glas- 
gow method , by ameliorating the overlap prob- 
lem, and to the imaginary chemical potential!^, 
which does not systematically bias the ensemble 
and thus needs large response to infinitesimal /i 
changes near ^ = 0. 

The fiuctuations of baryons are measured by 
the baryon number susceptibility 

X(T, fi) = dpifi, T)/9/i = d'logZifi, T)/ ^ 

Lattice results [|j indicate the range where 
x(T, /X = 0) is significantly different from zero : 
this is the candidate region for performing finite 
density calculations below Tc - a narrow, but not 
minuscule interval. 

As Z{iJ,) (we shall explicit only the /i depen- 
dence) is an even function of /i, we have: 

logZiij) = + a/i2 + + (2) 
logZ{iy) = K - av^ +b{v'^)+0{v^) (3) 

where v = i/i. Note that a — x{T, 0). Consider 

ae//(M) - ilogZi^i)-logZ{l^))/2^,^ (4) 
beffifi) = ilogZ{^i) + logZ{i^)^2K)/2fi^ (5) 

A fourth order polynomial describes the data 
where Oe// and beff do not depend on /i. 

The strong coupling limit of QCD displays con- 
finement and chiral symmetry breaking, and, be- 
sides having produced interesting, semiqualita- 



tive results on the phase diagram, offers a nice 
test bed for numerical methods. 

Consider the effective potential Vef / (related to 
Z a.s Z — {f VeffdX)^^ ) in the q = oo limit of 
QCD @: 

VeffiX,fi) = 2cosh{rNtNctJ.) + (6) 
sinh[{Nt + l)NcX\/sinh{NtX) 

where the variational parameter A (essentially, 

< ipip >) is to be determined by a minimum 
condition, Nt is the number of points in time di- 
rection, Nc the number of colors (which we fix: 
Nc ~ 3), r the asymmetry factor Ur/ap, and rNt 
is the inverse temperature. Note that formally 
T4//(A,/i) is an one dimensional partition func- 
tion which "remembers" its four dimensional ori- 
gin via A. 

For a purely imaginary chemical potential, 
cosh{rNtNc^J,) cos{rNtNcij)- We see the ex- 
pected periodicity 2tt /{rNtNc), and we note that 
the chemical potential term can be ignored for 
large Nt (zero temperature). Indeed, we have 
verified that in this limit any dependence on the 
chemical potential is lost jij, Fig.l : at T = 

< > as a function of (complex) /i is a con- 
stant in the half plane 3?(/i) < /ic, being the 
(real) critical chemical potential. By increasing 
the temperature, the effective potential changes 
with imaginary /x, and the chiral condensate in- 
creases with imaginary chemical potential, as it 
should. This behaviour is shown in Fig. 2, where 
we plot the chiral condensate as a function of a 
real, and a purely imaginary rather close to 
Tc- We note that the fourth order polynomial be- 
haviour anticipated above holds in a rather large 
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Figure 1. < ijiil) > as a function of real and imag- 
inary chemical potential (a constant) in the zero 
temperature limit. The only relevant quantity - 
pLc- is not amenable to an analytic continuation 
from imaginary /i. 



Figure 2. < > as a function of real (crosses) 
and imaginary chemical potential (plus'es) for T 
slightly below T^.. A polynomial representation 
(eqs. (2) and (3)) describes well either data set. 
a and b are computed as in eqs. (4) and (5). 



range of chemical potential, and that a simple 
second order polinomial works fine as well. 

In the same hot region reweighting too might 
prove useful. To check that, we have tested the 
Glasgow method in one dimensional QCD[||, a 
solvable model without SSB, but with baryons, 
whose partition function is formally identical to 
the effective potential above. The crucial differ- 
ence comes from the meaning of the parameter A 
which here is simply sh~^m. The bare mass m 
plays the role of an explicit breaking term, and 
< li^il) when m 7^ 0. Fixing iVc = 3, the rel- 
evant parameter is iV^A, which is the analogous 
of an inverse temperature. 

At large "temperature" (small NtX) we have 
verified that the Glasgow method reproduces the 
exact results 0, thus supporting the idea that 
reweighting methods can be successfully used in 
that regime. We would like to stress that this is 
a non trivial results since at large temperature 
quenched and exact results are distinct This 
contrasts with the low temperature regime of the 
same model, where the quenched approximation 
is apparently exact, and the success of the Glas- 



gow method is then guaranteed a priori (as the 
Glasgow method reproduces the quenched results 
for small statistics). 

Consider the zeros of the determinant, which 
control the singularities of the quenched model, 
and those of the partition function, the Lee 
Yang zeros controlling the singularities of the full 
model, in the complex plane e^. These Nt x 
zeros can be easily calculated using the results 
of and are shown in Figs. 3. The Glasgow 
method (Fig. 4) reproduces them correctly, as 
expected given that Z itself approaches the cor- 
rect value as mentioned above. The smoothening 
of the transition in the full model noted in ||] can 
be clearly read off Figs. 3, 4 : the zeros of Z move 
away from the real axis. It should be noticed, 
however, that in this case there is no pathological 
onset (no Goldstone particle here!): the Glasgow 
method only needs rearranging zeros on circular 
patterns. It is reasonable to hope that the same 
happens in high temperature, four dimensional 
QCD. From Fig. 3 we also understand why the 
quenched approximation is apparently working 
on long lattices 0]: the zeros will become dense 



3 




Figure 3. Distribution of the zeros of the determi- 
nant (100 configurations) and the exact results for 

the zeros of the partition function Z with Nf = A 
(diamonds) for one dimensional QCD. 



Figure 4. The zeros of Z =< Det > from the 
Glasgow method approach the exact pattern Fig. 
3. Correspondingly, number density and chiral 
condensate become nearly exact. 



on the same circle where the zeros of the deter- 
minant are randomly distributed. An important 
difference between the full and quenched transi- 
tion however remains: only the zeros of the par- 
tition functions are fixed by symmetries. Hence, 
only the nearest zero of Z (and not that of the 
determinant) has a calculable distance from the 
real axes ~ iJ.c/{3Nt), fulfilling a Lee- Yang-like 
scaling towards the Nt = oo singularity. 

We have given examples of clear signals ob- 
tained by use of an imaginary chemical poten- 
tial, and successful calculations with the Glasgow 
method in toy models related with QCD. We have 
suggested a possible simple line of analysis which 
combines and cross checks results from both ap- 
proaches. 
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